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The tidal interaction of a (rotating or nonrotating) black hole with nearby bodies produces changes 
in its mass, angular momentum, and surface area. Similarly, tidal forces acting on a Newtonian, 
viscous body do work on the body, change its angular momentum, and part of the transferred 
gravitational energy is dissipated into heat. The equations that describe the rate of change of 
the black-hole mass, angular momentum, and surface area as a result of the tidal interaction are 
compared with the equations that describe how the tidal forces do work, torque, and produce heat 
in the Newtonian body. The equations are strikingly similar, and unexpectedly, the correspondence 
between the Newtonian-body and black-hole results is revealed to hold in near-quantitative detail. 
The correspondence involves the combination /c2T of "Love quantities" that incorporate the details 
of the body's internal structure; ki is the tidal Love number, and r is the viscosity-produced delay 
between the action of the tidal forces and the body's reaction. The combination for is of order 
GM/c 3 for a black hole of mass M; it does not vanish, in spite of the fact that &2 is known to vanish 
individually for a nonrotating black hole. 



I. INTRODUCTION AND SUMMARY 
This work and its context 

Flanagan and Hinderer have recently shown [TJ |5] that 
the tidal interaction of two neutron stars in the inspi- 
ral phase of their orbital evolution leads to measurable 
effects in the gravitational waves emitted by the binary 
system. Such a measurement can be used to constrain 
the radius of each neutron star, and can thus reveal in- 
formation about the equation of state of nuclear matter. 
This exciting prospect has provoked a flurry of activity 
that aims to advance our understanding of tidal interac- 
tions between strongly self-gravitating bodies. 

Damour and Nagar [3J, and independently Binnington 
and Poisson [3], have formulated a relativistic theory of 
tidal Love numbers that applies to a neutron star (or 
any other type of compact object) that is slightly de- 
formed (from a spherical shape) by tidal forces. A tidal 
Love number is a dimcnsionlcss constant of proportion- 
ality that appears in the relationship between the tidal 
forces and the induced body deformation, measured by 
the mutipole moments of the mass distribution. Two 
types of tidal Love numbers occur in general relativity: 
An electric-type Love number associated with gravito- 
electric tidal forces, and a magnetic-type Love number 
associated with gravito-magnetic tidal forces. In each 
case the Love number refers to the shared multipole or- 
der of the tidal field and the induced deformation; there 
are two Love numbers for each multipole order. 

The electric-type Love number is directly analogous to 
the tidal Love number of Newtonian gravity (see, for ex- 
ample, Ref. [5]), and indeed, its relativistic definition re- 
duces to the Newtonian definition when the compactness 
C := 2GM/(c 2 a) of the body is small. (Here M denotes 
the body's mass, and a is its unperturbed radius; the 
compactness parameter C < 1 measures the importance 
of relativistic effects on the body's internal structure.) 



The magnetic- type Love number, on the other hand, is a 
relativistic quantity that vanishes in the Newtonian limit 
C — > 0. The Binnington-Poisson theory [1] applies just 
as well to black holes, and these authors have shown that 
the relativistic tidal Love numbers of a nonrotating black 
hole are all zero. 



Other notions of Love numbers have been introduced in 
Newtonian gravity [5] , and these also can be promoted to 
a relativistic framework. For example, Damour and Na- 
gar [3J have introduced "shape Love numbers" in the rel- 
ativistic theory of tidally deformed neutron stars; these 
relate the tidal forces to the displacement of the stellar 
surface. In a subsequent publication, Damour and Le- 
cian [5] extended the definition of shape Love numbers to 
black holes, generalizing a previous analysis by Fang and 
Lovelace [7]. These works have contributed significantly 
to our understanding of relativistic tidal interactions. 



My purpose in this paper is to contribute even further 
to this understanding. I show that certain aspects of 
the tidal interaction of black holes can be interpreted in 
purely Newtonian terms by exploiting a beautiful anal- 
ogy between the relativistic physics of black holes and 
the Newtonian physics of viscous fluids. This analogy was 
first noticed by Hartle [SJ [§] , and it was more fully fleshed 
out in the book Black Holes: The Membrane Paradigm 
|10j . Here the analogy is pushed even further, and un- 
expectedly, it is revealed to hold in near-quantitative de- 
tail. The analogy involves an additional "Love quantity" 
beyond the tidal and shape Love numbers described pre- 
viously. The new quantity r is not dimcnsionless; it pos- 
sesses the dimension of time, and in Newtonian physics it 
corresponds to the viscosity-produced delay between the 
action of the tidal forces and the body's reaction. I shall 
refer to r as the "viscous delay." 
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Flux formulae for black holes 

The physics to be interpreted concerns the rates at 
which a black hole changes its mass, angular momen- 
tum, and surface area as a result of the tidal interaction 
with external bodies. These were derived in Ref [TT1 112) . 
and the flux formulae are expressed in terms of quanti- 
ties £ ab that provide a characterization of the hole's tidal 
environment. This Cartesian 3-tensor is symmetric and 
tracefree (STF), and it represents the quadrupole mo- 
ment of the tidal field acting on the black hole. (Here 
and below the multipole expansion of the tidal field is 
truncated to the leading, quadrupole order; generaliza- 
tion to higher multipole orders is possible.) 

In Newtonian theory (see, for example, Ref. [13] for 
a modern treatment), the tidal- moment tensor is de- 
fined by first decomposing the Newtonian potential as 
U = £/body + C4xt, in which U\,o&y is associated with the 
reference body and U ex % is produced by the external ob- 
jects. We erect a Cartesian coordinate system x a that is 
attached to the center-of-mass of the reference body, and 
we define £ ab (t) := —d ab U cx t, in which the spatial deriva- 
tives of the external potential are evaluated at x a = 0, 
the position of the center-of-mass. Because Z7 ex t satis- 
fies Laplace's equation in the body's neighborhood, £ ab 
is tracefree in addition to being symmetric in its indices. 

We consider a context in which the tidal forces ex- 
erted by the external bodies are weak, and in this case 
it is appropriate to express the external potential as 
U cxt = U (t) + g a {t)x a - \£ ab {t)x a x b + • ■ • , a Taylor ex- 
pansion about the body's center-of-mass. The spatially- 
constant term Uo(t) is irrelevant, and g a (t) :— d a U cx t is 
responsible for the body's acceleration in the field of the 
external objects; the remaining terms are responsible for 
the tidal forces acting on the body. 

In general relativity, the tidal moments £ a b{v) '■= CoaOb 
are defined in terms of the components CoaOb of the space- 
time Weyl tensor evaluated in a region of spacetime that 
is far away from the black hole, but still close relative 
to the external bodies. In this case the indices refer to 
an inertial frame moving with the black hole, and v is 
an advanced-time coordinate on the event horizon. (See 
Refs. PUCE2] for details.) 

For a nonrotating black hole of mass M, the rates at 
which the tidal interaction changes the hole's mass and 
angular momentum are given by Eqs. (8.38) and (8.39) 
of Ref. fTT]; these are 



M = ^M e £ ab £ ab 
45 



and 



j=-%M\e» cd £ 



-b„d- 



£ab 



(1.1) 



(1.2) 



respectively. Here an overdot indicates differentiation 
with respect to advanced-time v, e a bc is the Cartesian 
permutation symbol, and J := J a s a is the component 
of the vector J a in the direction of the (arbitrary) unit 



vector s a . The rate of change of the surface area is 
obtained from the first law of black-hole mechanics ap- 
plied to a nonrotating black hole, (k/8it)A — M, where 
k = (4Af) _1 is the hole's surface gravity. I employ rela- 
tivistic units, so that G = c = 1. 

For a rapidly rotating black hole of mass M and 
angular-momentum J, the rate at which the angular mo- 
mentum changes as a result of the tidal interaction is 
given by Eq. (9.39) of Ref. [II]; this is 



J = 



2(l+3 X 2 ){£ab£ ab 



Z(l+ l l x *)(£ ab3 »£ 



bea s c\ 



15 



+ ^ x l(£ abS a s bf 



(1.3) 



Here \ := J/M 2 is the dimensionlcss Kerr parameter 
(which is limited to the interval < X < 1)> the unit 
vector s a is the direction of the hole's rotation axis, and 



O 



X 



H 



2M(1 + y/l - x 2 ) 



(1.4) 



is the angular velocity of the event horizon. An equa- 
tion for M is not available in this case, but at the level 
of accuracy achieved by Eq. (1.3), A can still be ob- 
tained from the first law: (k/8tt)A = —VLhJ, where 
k = (2M) _1 yl — x 2 (l — \/l — X 2 ) 1 is the surface grav- 
ity of a Kerr black hole. 

Equations (1.1 1, (1.2|, and (1.3 1 differ from the orig- 



inal expressions displayed in Ref. |llj by terms involv- 
ing the gravito-magnetic tidal moments B ab {v). These 
are switched off because they possess no Newtonian ana- 
logues and cannot, therefore, be involved in a correspon- 
dence with the tidal dynamics of a Newtonian body. The 
neglect of the gravito-magnetic terms can be motivated 
on the basis of a post-Newtonian approximation of the 
tidal fields [13] ; in this context the terms in M and J that 
involve B ab are suppressed relative to those involving £ ab 
by factors of order (iWb/c) 2 <C 1, where v or b measures 
the orbital velocity of the external bodies. 



Tidal coupling of a Newtonian, viscous body 



The key to a Newtonian interpretation of Eqs. (1.1 1 



(1.2), and (1.3 1 is the relation (Sec. Ill D) 



Qjk 



° 5 ( 



r£ jk 



(1.5) 



that holds between the mass quadrupole moment Qjfc of 
a Newtonian body of mass M , unperturbed radius a, and 
angular velocity fin, and the quadrupole moment £jk of 
the applied tidal field (an overdot indicates differentiation 
with respect to time t). This relation is formulated in a 
frame x J that rotates uniformly with angular velocity 
relative to the global inertial frame (the frame corotates 



3 



with the body); this differs from the original frame x a , 
whic h mo ves with the body but does not rotate. Equa- 



tion (1.5) involves the STF tensor Cjk '■— ^Sjk ~ SjSki 
where Sj is a unit vector that points in the direction of the 
rotation axis. It involves also the dimensionless numbers 
n 2 (rotational Love number) and k 2 (tidal Love number), 
and the viscous delay r. The body's mass quadrupole 



moment is defined by Q 



J p(x 3 x k 



2 <p' fe ) (Pi 



which p is the mass density; this STF tensor vanishes 
when the density profile is spherically symmetric. 

The first term on the right-hand side of Eq. ( |1.5[ ) de- 
scribes the body's rotational deformation (flattening of 
the poles, bulging of the equator). The scalings with Vt 2 H 
and a 5 are dictated by dimensional analysis, and the ro- 
tational Love number n 2 incorporates the details of the 
body's internal structure. The second and third terms 
collectively describe the tidal deformation. Once more 
the scaling with a 5 is dictated by dimensional analysis, 
and the tidal Love number k 2 depends on the body's in- 
ternal structure. The third term is contributed by the 
body's internal viscosity, and dimensional analysis re- 
veals that r oc au/M, with v denoting the (averaged) 
kinematic viscosity of the fluid. 



The relation of Eq. (1.51 relies on an assumption that 



the time scale associated with changes in £jk(t) is long 
compared with the viscous delay r. In this context the 
second and third terms can be packaged approximately 
as — ^k 2 a 5 £jk(t — r), which reveals that viscosity does 
indeed introduce a delay between the action of the tidal 
forces (described by £jk) and the body's reaction (mea- 
sured by Qjk)- The viscous delay produces a misalign- 
ment between the figure of the body deformation (de- 
scribed by Qjk) and the figure of the tidal forces (mea- 
sured by £jk)- In the case in which the tides are produced 
by a single external body, the tidal bulge points in a di- 
rection that is not quite aligned with the direction of the 
external body. 



Equation (1.51 holds in the body's rotating frame x 3 



and Ejk is the rate of change of the tidal moment as 
perceived by an observer corotating with the body. In 
the nonrotating frame x a the quadrupole moment is given 
instead by (Sec. Ill E) 



2 2 - 

Qab = -n 2 a b tf H C ab --k 2 a b [£ ab -T£ ab -TA£ ab ). (1.6) 



In this expression C a b '■= \5 ab — s a s b , in which the unit 
vector s a continues to point in the direction of the ro- 
tation axis; £ ab is now the rate of change of the tidal 
moment as perceived by a nonrotating observer, and the 
last term 



A£ ab := 2^l H e cd ^ a £ c b ^s d 



(1.7) 



accounts for the rotation of the body frame relative to 
the global inertial frame. 



Flux formulae for Newtonian body 

The rate at which the tidal forces do work on the New- 
tonian body is calculate d as W — \Q ab £ ab (Sec. IV A), 
and substitution of Eq. (1.6 1 yields 



W= l -{k 2 T)aH ab (£ ab + A£ ab ) 



(1.8) 



Notice that W is proportional to k 2 r: No (net) work is 
done unless the body is deformed by the tidal forces, and 
no (net) work is done unless the deformation is delayed 
with respect to the application of the forces. In the ab- 
sence of viscosity, the tidal bulge points in the direction 
of the external body, and the transfer of gravitational 
energy from the tidal field to the body is fully reversible 
and does not lead to a net performance of work. In the 
presence of viscosity, on the other hand, the misalign- 
ment of the tidal bulge allows for an irreversible transfer 
of energy from the tidal field to the body, as described 
by Eq. Oh. 



The torque exerted on the body by the tidal forces 
produces a change of angular momentum describe d by 
j a = - £abc Qb p £pc (Sec. IV B). Substitution of Eq. Q | 
and projection along s a yields 

n H J = ~(k 2 T)a">A£ ab (£ ab + A£ ab ). (1.9) 

The result is multiplied by f2# to convert J into E TOt = 
Q,H J, the rate of change of the body's rotational energy. 
Once more the result is proportional to k 2 r: The angular 
momentum changes if and only if the body is deformed 
by the tidal forces and the deformation is delayed with 
respect to the application of the forces. In the absence 
of viscosity, the tidal bulge points in the direction of the 
external body, and the tidal forces have no opportunity 
to exert a torque. In the presence of viscosity, on the 
other hand, the misalignment of the tidal bulge allows 
the tidal forces to exert a torque, and the end result is 
described by Eq. (1.9 1. 



The rate at which viscosity generates heat can be de- 
termined from the energy-balance e quat ion ( Sec. IV D) 
Q = W-Cl H j. Substitution of Eqs. Ol and Esfy yields 



Q=-(k 2 T)a 5 (£ ab + A£ ab )(£ ab 



A£ ab 



(1.10) 



Newtonian interpretation of black-hole fluxes 

The black-hole and Newtonian-body results can be 
put in a close corresponden ce if we associate the W of 
Eq. (L8| with the M of Eq. a nd the J of Eq. |l~9 | 

with the J of Eqs. (1.2 1 and (1.3 1. With these associ- 
ations, the energy-balance equation Q = W — QhJ be- 
comes the first law of black-hole mechanics, (k/8tt)A = 
M — HhJ, when we also associate Q with (k/87t)^4; this 
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association is natural if we think of k/2it as the black- hole 
temperature and \A as its entropy. 

The nonrotating black hole corresponds to the case in 
which £ ab dominates o yer A£ n h in the Newtonian equa- 
tions; recall from Eq. (1.7 1 that A£ ab is proportional to 



£Ih and therefore vanishes when the body is nonrotating. 
In this case Eqs. (1.8 1 and (1.9 1 give 



1 



M = -(k 2 r)a 5 £ ab £ ab , J = --{k 2 T)a 5 (e a cd £ cb s d )£ ab , 

(1-11) 

and these equations bear a striking resemblance with 
Eqs. (O) and pi. 



The rapidly rotating black hole corresponds to the 
case in which A£ ab dominates over £ ab ; the changes 
in the tidal moment are produced almost entirely by 
the body's own rotation. In this case Eq. (1.9 1 gives 



Cljfj = — ^(k 2 r) a 5 A£ ab A£ a i,, and expanding the tenso- 
rial expression produces 



2(£ ab £ ab ) - 3(£ ab s b £ a c , 



(1.12) 



Comparison with Eq. (1.3) reveals another striking sim- 
ilarity, especially if we neglect the terms of order \ 2 
in E q. (|1.3 |. This is appropriate, because according to 
Eq. \l.ty , x 2 c an be thought of as a relativistic correction 
of order (u rot /c) 2 , in which v TO t = afln stands for the ve- 
locity of a fluid element on the surface of the rotating 
body. 

Closer examination reveals that the correspondence 
produces agreement between all numerical coefficients 
when we make the assignment 



(k 2 T) 



5 _ 16 / GM\ / GM\ 



15 V c 3 J\ c 2 J 



(1.13) 



in the Newtonian equations. (I have restored the factors 
of G and c, which were previously set equal to unity). 
The equation indicates that as might be expected, the 
horizon length scale GM/c 2 must play the role of the 
body radius a, while the horizon time scale GM/c 3 must 
play the role of the effective viscous delay k 2 r. This 
last result implies that the horizon can be assigned an 
effective kinematic viscosity k 2 v ~ GM/c. 

The picture that emerges from this comparison is one 
in which the event horizon behaves as if it were a ficti- 
tious membrane of viscous fluid; the fluid elements move 
with the speed of light, and the streamlines are identi- 
fied with the horizon's generators. This view, of course, 
is compatible with the membrane paradigm of black-hole 
physics [TU] , 

The correspondence between the black-hole and 
Newtonian-body fluxes depend on the combination k 2 r 
of "Love quantities," which is constrained by Eq. (1.131; 



the constraint applies to rapidly-rotating as well as non- 
rotating black holes, and it implies that k 2 r must be 
a finite (nonvanishing) quantity. It is known, however, 
that k 2 = for a nonrotating black hole [1]. One must 



therefore think of the combination k 2 r as being insep- 
arable for (nonrotating) black holes. While k 2 vanishes 
individually, k 2 r is nevertheless a finite quantity of order 
GM/c 3 . It is not helpful to think of t as an individuated 
quantity that happens to be infinite. 



Organization of the paper 

In the remaining sections of the paper I elaborate the 
Newtonian theory of the interaction of a fluid body with a 
tidal environment characterized by the tidal quadrupole 
moment £ ab (t). In addition to establishing other inter- 
esting results, I provide derivations for Eqs. (1.5 1, (1.6 1, 



(1.8), (1.91, and ( 1.10 1 



I begin in Sec. II with a warmup exercise in which the 
fluid's viscosity is set equal to zero; viscosity is intro- 
duced next in Sec. III. In Sec. IV I calculate the rates 
at which the tidal forces do work on the body, change 
its angular momentum, and lead to the production of 
heat. For simplicity I perform all calculations assuming 
that the body consists of an incompressible fluid; general 
results are obtained from the special cases by inserting 
Love numbers where appropriate. 



II. TIDAL INTERACTION OF A NEWTONIAN 
BODY; NONVISCOUS CASE 

We consider a body of mass M, radius a, and angu- 
lar velocity fin that would be spherical in isolation and 
in the absence of rotation. The body is made of an in- 
compressible fluid, it rotates rigidly, and it is put in the 
presence of external objects that exert tidal forces. As 
a result of its rotation and tidal interaction, the body 
acquires a deformation. We calculate this deformation, 
first assuming that the fluid's viscosity can be neglected. 



A. Scales and assumptions 

We already introduced the scaling quantities M, a, and 
Qh- The tidal interaction is characterized by an external 
mass scale M' and a distance scale b; it is measured by 
the tidal moment tensor £jk, whose scale is given by £ ~ 
GM' /b 3 . The time scale associated with changes in the 
tidal field, as viewed in the global incrtial frame, is 
where f2 is an angular velocity associated with the motion 
of the external objects. In the body's rotating frame the 
time scale becomes |f2 — f2jj| _1 , and the tidal field is 
static if the angular velocities match at all times. The 
tidal forces produce a displacement of a fluid clement at 
the body's surface; the length scale of this displacement 
is £ ~ a 4 £/(GM), and we assume that (<a. We also 
assume that the deformation produced by the rotation is 
small, but we make no assumption regarding the relative 
size of these effects. 
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The velocity scale associated with fluid motions in the 
rotating frame is v ~ £|fi — f2jy| _1 , or c ~ a 4 |f2 — 
tin \~ 1 £/(GM ). We assume that v is much smaller than 
the linear velocity af2# of a fluid clement following the 
body's rotation (as viewed in the global incrtial frame). 
The inequality v <C aftu implies that the Coriolis ef- 
fect can be ignored, but centrifugal terms will make an 
appearance in the fluid equations. This assumption con- 
strains the size of the tidal interaction. Assuming that 
Cljj is larger than or at least comparable to CI, the in- 
equality translates to £ < GM/a 3 , or (a/b) 3 < M/M'; 
the body must be well separated from the objects that 
produce the tidal forces. As an additional assumption 
we take v\Q — fin the velocity scale divided by the 
tidal time scale, to be much smaller than £a 7 the scale 
associated with the tidal forces (per unit mass) acting 
within the body. This inequality allows us to neglect 
the inertial term dtv in the fluid equation. It translates 
to (f2 — Qh) 2 <C GM/a 3 , and the assumption constrains 
the size of the angular velocities; they must be small com- 
pared with the Keplerian angular velocity of a particle at 
the body's surface. All the assumptions made here are 
reasonable, and they do not severely restrict the range of 
physical situations that can be considered. 



B. Fluid equations; nonrotating frame 

We first formulate the fluid equations in a frame x a 
that moves with the body's center-of-mass but does not 
rotate relative to the global inertial frame. (This frame 
was denoted x a in Sec. I, and here we put an overbar on 
the index to increase notational clarity.) The frame x a 
is translated by r a {t) relative to the origin of the global 
inertial frame; this is the position of the body's center- 
of-mass. 

The equations of fluid dynamics are reviewed, for ex- 
ample, in Sees. 1 and 2 of Ref. |15j . A perfect fluid is 
governed by (i) the continuity equation 

d tP +d- a (pv a ) = 0, (2.1) 

in which p is the mass density and v a the velocity field; 
(ii) Euler's equation 

dv a 

p—7- = -d a p + pdaU - prs., (2.2) 
at 

in which dv a /dt = d t v a + v b d} ) v a is the advective time 
derivative, p the pressure, U the gravitational potential, 
and the last term is an inertial force density associated 
with the translation of the frame x a relative to the global 
inertial frame; (iii) Poisson's equation 

V 2 U = -AirGp (2.3) 

for the gravitational potential; and (iv) an equation of 
state that relates the pressure to the density. 



In the case of an incompressible fluid, p = constant 
and the fluid equations become 

d- a v a = (2.4) 

and 

d ^ = d- a (-p/p+U-r l x l ). (2.5) 

In this case the pressure is unrelated to the density. 

It is useful to decompose the Newtonian potential as 
U = [/body + U ex t, in which [/body is associated with the 
body and C/ ex t is produced by the external objects. We 
also decompose the external potential as U cxt (t,x a ) — 
U (t)+ga(t)x a +U tidal {t,x a ), in which U (t) := U cxt (t,r a ) 
and g a '■— <9 a C/ cx t(i, r a ); this provides a definition for the 
tidal potential J/tidai- Because the body moves according 
to f s = g a , the final form of Euler's equation is 

dv & 

= 9 a {-P/P + ^body + ftidal) • ( 2 -6) 

C. Fluid equations; rotating frame 

The transformation from the frame x a = [x, y, z] to the 
rotating frame x^ = [x, y, z] is described by 

x — xcosflnt + ysmflift, (2.7a) 
y = —x sin ilnt + y cos tint, (2.7b) 
z = z, (2.7c) 

if we take the rotation axis to coincide with the z-axis. 
The transformation can be written more formally as 

A '„..<"••. (2.8) 

The velocity vector transforms as 

v x = v x cos Vt H t + v° sin £l H t 

+ xsintlnt + ycosilnt), (2.9a) 
v y = —v x sin flfft + v ^ cos t 

— flH(x cosflnt + ysmfl H t^, (2.9b) 

v z = V s , (2.9c) 

which can be written formally as 

v 3 = A\(v a + Av a ), (2.10) 

where 

Av a :=T a - b x h (2.11) 

accounts for the time-dependence of the transformation. 
The matrix T a - b is formally defined by T a - b = A^ A^, 

where A"^ is the inverse of the matrix A 3 a . Calculation 
reveals that 

n = n H e%-S 5 , (2.12) 
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where e a j g is the permutation symbol and the vector s c — 
[0,0, 1] denotes the direction of the rotation axis. 

Making the substitution in Eulcr's equation produces 
additional terms that originate from the nonincrtial na- 
ture of the transformation. The first collection of ad- 
ditional terms are linear in both flu and v J , and are 
responsible for the Coriolis effect; these we neglect, ac- 
cording to the assumptions formulated in Sec. II A. The 
second collection of terms are quadratic in and inde- 
pendent of , and give rise to centrifugal repulsion; these 
we keep. To simplify the equation further we neglect all 
terms quadratic in v° and the inertial term dtv j , again 
in accordance with the stated assumptions. 

As a result of these substitutions we find that the fluid 
equations in the rotating frame are djv^ = 0, 

dj (-p/p + U body + C/ tida i + fl 2 H C) = 0, (2.13) 

and the gravitational potentials are determined by 
V 2 [/ b ody = -47rGp and V 2 [/ t idai = 0. Here Q 2 H C = 
\Q 2 H {x 2 + y 2 ) is the centrifugal potential. It is useful to 
express C as 

C= l -r 2 + l -C jk x^x k 1 (2.14) 

where 

C jk = ^5 jk - SjSk (2.15) 

is a symmetric-tracefree (STF) tensor. The first term in 
C's decomposition is spherically symmetric; the second 
has a quadrupolar (£ = 2) structure. 

D. Tidal field 

The gravitational potential was decomposed as 

U = £/b od y + Utidal, (2-16) 

with [/body representing the potential of the (deformed) 
body, and [/tidal representing the tidal potential. To lead- 
ing order in a Taylor expansion of the external potential 
about the body's center-of-mass, this is 

U tida . l = ~£ jk {t)x j x k , (2.17) 

where the tidal moment £j k is a time-dependent STF 
tensor that does not depend on spatial position. The 
tidal potential also has a quadrupolar structure. 

A specific example of a tidal field is one produced by an 
external object of mass M' on a circular orbit of radius b 
around the reference body. If we let m? be a unit vector 
that points in the direction of the external object (from 
the body's center-of-mass) , a simple computation reveals 
that the tidal moment is given by 

GM' / 1 \ 

£ jk{t) = — gp- \m m k - ^S jk j. (2.18) 



If we take the external object to move in the equatorial 
plane of the rotating body, the unit vector is given by 

m j (t) = [cos(n - n H )t,sm(n - n H )t,o] (2.19) 

in the rotating frame, where ft is the orbital angular ve- 
locity. 

E. Body deformation 

The tidal forces and body rotation produce a deforma- 
tion in the body's shape, which can be described by 

R = a(l + e jk n j n k ), (2.20) 

where -n? — x^ jr is the radial unit vector. The equation 
r = i?(n J ) describes the position of the body's deformed 
surface. We let the deformation assume a quadrupolar 
shape, in accordance with the fact that each deforming 
force is quadrupolar in nature; ej k is a STF tensor that 
must be determined by solving the fluid equations. 

F. Body potential 

The body potential is given by 

[/ - GM + Vo x ° xk 

<^body,out — 1 ^LrVj'fe ~ (Z.ZL) 

outside the body, and 

[/body.in = U (r) + U jk x j x k (2.22) 

inside the body. In each case the deviation from a spher- 
ical potential is measured by a STF tensor; we have the 
quadrupole-moment tensor Qj k outside, and Uj k inside. 
The body mass is related to the density by 

Air 

M = —pa 3 . (2.23) 

The quadrupole-moment tensor is determined by evalu- 
ating Q jk = J px( jk ^ d 3 x, where := x j x k - \r 2 5 jk \ 
a quick computation returns 

Q jk = 2 -Ma 2 e ]k . (2.24) 

The internal potential Uq is obtained by solving Poisson's 
equation for a spherical configuration of constant density 
p and radius a; the result is 

£MrH^(3-r 2 /a 2 ). (2.25) 

And finally, Uj k is obtained by demanding continuity 
of the body potential across the deformed surface; this 
yields 

3 GM 

l - iri (2-26) 
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G. Pressure 



The pressure field inside the body is expressed 



as 



p/p = Po(r)/p + Pjk 



X^ x^ 



(2.27) 



where pjk is another STF tensor. This form is mo- 
tivated by the fact that according to Eq. ( [2.131 ), the 
pressure must satisfy the Poisson equation W z (p/ p) — 
~4irGp + VL 2 H . Because the right-hand side is spherically 
symmetric, po/p must satisfy this equation by itself, and 
the additional term must be a solution to Laplace's equa- 
tion; a quadrupolar structure necessarily selects the as- 
sumed form pjkX : 'x k . 

The spherically-symmetric component of the pressure 
is easily obtained by integrating Poisson's equation; we 
get 



Po(r)/p 



GM 
2a 



(2.28) 



when we demand that the unperturbed pressure vanish 
on the unperturbed surface r = a. An expression for pjk 
is obtained when we demand that the perturbed pressure 
vanish on the perturbed surface. If we let £ := R — a = 
aejkn : 'n k describe the deformation, we have that = 
p(a + £) = po(a) + ap' (a)eji,n : >n k + a 2 pj^-n?n k , in which 
a prime indicates differentiation with respect to r. The 
computation returns 



GM 



Pjk 



&jk 



(2.29) 



after neglecting the correction of order VL 2 H /(AiiGp) in 
Po{r). 



H. Solution 



Euler's equation (2.13) produces the constraint 

- Pjk + U jk - \£j k + \n 2 H C 3k = (2.30) 

when we insert Eqs. Q2.14[ ), \2.\7\ , p.22) , and \2.27\ into 
it. Since pjk and Ujk are both proportional to the defor- 
mation tensor ejk, we quickly deduce that 



5 a 6 



ejk 



\fl 2 H Cjk - £jk)- 



(2.31) 



4 GM 

This result informs us that the deformation is a combi- 
nation of rotational and tidal effects, as expected. We 
also have 

p jk = -{n 2 H Cjk - £jk) 



u 



(QjrC 



H^jk 



£ 



and the quadrupole-moment tensor is 
1 a 5 



Ijk 



2 G 



(n 2 H G 



£ 



(2.32) 
(2.33) 

(2.34) 



III. TIDAL INTERACTION OF A NEWTONIAN 
BODY; VISCOUS CASE 

We now incorporate viscosity into the fluid model. The 
fluid's shear viscosity is measured by v, the coefficient 
of kinematic viscosity, which we assume to be uniform 
within the fluid. 



A. Scales and assumptions 

We continue to deal with the scaling quantities intro- 
duced in Sec. II A, and with one exception, we make 
the same assumptions regarding them as we did in the 
preceding section. The only change concerns the time 
scale of the problem, which previously was provided by 
the orbital dynamics as viewed in the rotating frame; the 
only relevant time scale was |f2 — I -1 . The addition of 
viscosity introduces another relevant time scale, namely 
the delay between the action of the applied tidal force 
and the body's response. As we shall see, the time scale 
for the viscous delay is r ~ av/(GM), and we assume 
that this is much shorter than the orbital time scale: 
t \Q, — Oj^- 1 1 ; this inequality constrains the size of 
the kinematic viscosity. 

The change of time scale has repercussions on the ve- 
locity scale. It is now given by v ~ £/r, where £ ~ 
a 4 £ /(GM) is still the length scale associated with fluid 
displacements. To neglect the Coriolis effect with still de- 
mand that v <^ afln, but the inequality now translates 
to £ < (GM/a 3 )(fl H T), or (a/b) 3 < (M/M'){Q h t); 
since fijjr < 1 we find that the bodies must be even 
more widely separated. To neglect the inertial term in 
the fluid equation we now demand that v/t <C £a, and 
this translates tor> yj a 3 / (GM), another constraint on 
the kinematic viscosity. (The constraints are compatible 
with each other whenever \Q — f2#| <C t -1 <C Qk, where 
Q K = yjGM/a 3 is the Keplerian angular velocity of a 
particle at the body's surface.) 



B. Fluid equations 

Under the stated assumptions, the fluid equations in 
the rotating frame are now (see, for example, Sec. 15 of 
Ref. [H]) 



dji) 3 = 



(3.1) 



and 



dj (-p/p + U hody + C/ tidal + C) + jA7 V = 0, (3.2) 

with the last term giving rise to all effects associated with 
the viscosity. The centrifugal potential is still given by 



C - hl 2 H r 2 + Ul 2 H C jk x j x k , 



(3.3) 
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and the internal gravitational potential is still given by 

U hody + f/tidai = U (r) + U jk x j x k - ^£ jk x j x k , (3.4) 

with U (r) = GM(3 - r 2 /a 2 )/(2a) and U jk = 
3GMej k /(5a 3 ). The body deformation is still described 
by R = a(l + ej k n?n k ), and the quadrupole moment is 
still Q jk = |Ma 2 e jfe - 

The presence of viscosity affects the nature of the 
boundary conditions at the body's surface. Instead of 
the no-pressure condition imposed previously, we must 
now impose the no-normal-stress condition 



prij - pv(d J v k + d k Vj)n k = 



(3.5) 



at r = R. In addition, we must relate the normal com- 
ponent of the velocity vector to the rate of change of the 
surface displacement: 



d t R. 



(3.6) 



Unlike the preceding equation (which involves the un- 
perturbed piece of the pressure), this equation can be 
formulated at r = a instead of r = R, because it involves 
perturbed quantities only. 

Because the pressure still satisfies the Poisson equation 
V 2 (p//o) = — 4irGp + we may still express it as 



P/P = Po(r)/p + p jk x J x k . 



(3.7) 



The expression for po( r ) is unchanged, but pj k is no 
longer equal to GMej k /a 3 ; its new value will be deter- 
mined by the new boundary conditions. Using the avail- 
able information, we find that these become 



v(djV k + d k Vj)n = a n,j[p kn 
and 



GM 



e kn n n 



VjinP — aej k nPn k . 
Both equations are now formulated at r = a. 



(3.8) 



(3.9) 



C. Velocity field 

At this stage the unknowns are the velocity field w J , 
the displacement tensor ejk, and the pressure tensor 
Pj k ; these must all be determined in terms of Cj k and 
£j k . The velocity field obeys a Poisson equation with a 
quadrupolar source term. The solution must be a £ — 2 
vectorial harmonic, and it can therefore be expressed as 
a linear superposition of terms r 2 Vj k x k , V kn x^x k x n , and 
Vj k x k , with Vj k a time-dependent STF tensor; the super- 
position is constrained by the divergenc e-fre e condition 
djV J = 0. The boundary condition of Eq. (3.9 ) relates Vj k 



to ej k , and Eq. (3.8) determines pj k in terms of ej k and 
its time derivative. Finally, the Navier-Stokes equation 



(3.2) gives rise to a differential equation for the displace- 
ment tensor ej k . 

After going through these manipulations and sorting 
out the algebra, we arrive at the differential equation 



5 a" 



T€-j k ~r &jk — 



4 GM 



-{n 2 H G jk - Ejk), 



where 



19 av 



(3.10) 



(3.11) 



2 GM 

is the viscous delay. We also obtain the pressure tensor 



Pjk = 

and the velocity field 
1 



GM j 



4r 
95 



Gjk 



(3.12) 



5a 2 



(5r 2 -8a 2 )e jk x k -2e kn x J x k x n . (3.13) 



D. Solution 



The solution to the differential equation re + e = / is 



e(t) 



e(0) + - / e t,/T f(t')dt'}. (3.14) 

T JO J 



Two integrations by parts bring this to the form 
e(t) = [e(0)-/(0) + r/(0)] e - t ^ + /(i)-r/(t) 

+ re- t/T f e* ' /T f(t')dt> '. (3.15) 



At times ( » r we may neglect the first group of terms, 
and we may also neglect the integral if f(t) varies on a 
time scale that is long compared with r. Under these 
circumstances the solution is 

e(t) = f(t) - rf(t) + O[e(0)e-^] +0[r 2 f(t)}. (3.16) 

Within the stated error, this is e(t) = f(t — r). We see 
that the displacement at time t is related to the behavior 
of the driving force at the earlier time t — T, and the delay 
is given precisely by r. The viscosity therefore introduces 
a phase lag between the applied force and the response, 
and as we shall see, this dephasing turns out to be crucial 
for the physics of tidal work, torque, and heating. 

In our context the role of the driving force / is played 
by the tidal moment £j k , which varies on a time scale of 
order |f2 — We already have stated our assump- 

tion that this time scale is much longer than r. In these 
circumstances, and ignoring the transients that decay ex- 
ponentially, the solution to Eq. (3.101 is 



ejk 



5 a d 



4 GM 



■'jk 



■rSjk). 



(3.17) 
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With this we find that 
5 

Pjk = j 
U ih = - 

3 K A 



(p%Cjk — £jk + ' 



(^l^jCjk — £jk + T £jk 



(3.18) 
(3.19) 



and 



4GM 



(5r 2 - 8a 2 )£j k x k - 2£ kn x j x k x n ] . (3.20) 



The quadrupolc moment tensor is now 
2G 



Qjk = -^q [pH c jk - £jk + T~£jk 



(3.21) 



E. Transformation to the nonrotating frame 

For later purposes we transform our results from the 
rotating frame x 3 back to the nonrotating frame x a ; the 
transformation is described by x a = A a ^x 3 \ and it is the 
inverse of the one considered in Sec. II C. 

As we have seen, the velocity field transforms as 



where 



v 5 - = A a Jv j + Av j ), 



Av 3 := T\x k 



(3.22) 



(3.23) 



accounts for the time-dependence of the transformation; 
here 



yj — \] A" 
1 k ■— A a A k 



-Clue 



kn b i 



(3.24) 



with s™ = [0,0,1] denoting the direction of the rotation 
axis. 

In the nonrotating frame the tidal moment is given by 



£- al = A*- a A\£ jk 



(3.25) 



The transformation of its time derivative, however, must 
account for the time-dependence of the transformation. 
Here we have 

4 5 - A 3 - a A\(£ jk - A£ jk ) = A\A\ £ 3k - A£- ah (3.26) 
with 

A£jk := T n j£ n k + T n k £j n = 2U,H^ pq (j£ p k ^s q (3.27) 

and 

A£- al := A 3 a A k h A£ jk = 2Q H e- oS(a £ %/. (3.28) 
The transformation of the qu adrup olc-m omcn t tensor 



is Q- a i = A\A\Q jk . With Eqs. (3.21 1 and (3.261, this is 



Qab — {pH^ab &ab + r ^ab + T ^£ab j 



(3.29) 



with the last term accounting for the effect of the body's 
rotation on the tidally-induced quadrupole moment. 



F. Example 

To illustrate the way in which the preceding (sim- 
ple, but potentially confusing) tensorial transformations 
work, we examine the case in which the rotating- frame 



tidal moment is given by Eq. (2.181 
GM 



-'jk 



GM 1 / 1 \ , 

-gp- [mjirik - -S jk J , (3.30) 



with m? = [cos(f2 — f2#)t, sin(f2 — D,jj)t,0]. The trans- 
formed tidal moment is 



-{rriami - ^ 5 &)> (3-31) 



GM 
ah 3b 3 

with m a = [cosf2i,sinfM,0]. The time-derivative of the 
rotating-framc tidal field is 



£ 



jk 



GM' 
' 3b 3 



(fl - tt H )(mj<i> k + (pjm k ), (3.32) 



where (ft 3 :— [— sin(f2 — f2ff)i, cos(f2 — 0]. From this 

we find that 



A j -A k -£ i. 

Iy a 11 b 3 k 



GM' 
"363" 



(fl-Q H )(ma<j>b + <f>a,mi), (3.33) 



with <f) a — [— sin ilt, cos Sit, 0]. On the other hand, 
GM' 



A£ 



jk 



3b 3 



-Q H (mjfa + Mjirik) (3.34) 



and 



GM' 

A 4b = -jjp-fMmaifo + 0am-b). (3.35) 

We see that these pieces add up correctly to produce the 
expected 



£r,h 



GM' 
3b 3 



The combination of terms that appears in Q^, however, 
is 



ab 



GM' 
3b 3 



(fi - Q H )(ma<j)- b + <\>-ayn- h ) , (3.37) 



and we see that this vanishes when the external object is 
corotating with the body. 



IV. TIDAL WORK, TORQUE, AND HEATING 

A. Tidal work 

The rate at which the tidal forces do work on the ro- 
tating body is calculated as 

b\ „,a j3 - 



W= J {-p£- al x l )v°d 
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This is the integral of the tidal force density —p£ ab x b 
times v a , the rate of fluid displacement. This is equal to 



or 



(4.1) 



The total derivative corresponds to a change of a state 
function, and we shall ignore such (reversible) changes in 
this analysis. Our final expression for the rate of tidal 
work shall be 



(4.2) 



If we insert Eq. (3.29) in this we find that the £l H C &b 



and £ ab terms contribute total derivatives, which we con- 
tinue to ignore. The only nontrivial contributions come 
from the terms proportional to r; we get 



W=—£ ab {£- ab + A£- ab ). 



We recall that 



and 



A£- al = 2n H e 5d{ - a £% )S 



19 av 
~2GM' 



(4.3) 



(4.4) 



(4.5) 



The equation reveals that there is no net tidal work in 
the absence of viscosity. The dephasing between Q ab and 
£ ab created by the viscous delay is therefore an essential 
piece of the physics. In the absence of viscosity, the fluid 
displacement is always in phase with the applied force, 
and this leads to a reversible transfer of gravitational 
energy between the body and the tidal field; the work 
done necessarily averages to zero. 

For an external object of mass M' on a circular orbit 
of radius b and angular velocity fi, the rate of tidal work 



There is no work done if the external object is corotating 
with the body. 



B. Tidal torque 

The total torque exerted by the tidal forces is 



and this is equal to J a , the rate of change of the body's 
angular-momentum vector. Evaluation of the integral 
gives 



Ja — e abcQ p£ P 



(4.7) 



The angular-momentum vector can be decomposed as 
J a = Js a , in terms of a magnitude J and a direction s a . 
The change in angular momentum is then expressed as 
J a = Js a + Js a . The first term describes the change in 
the body's rate of rotation, while the second describes 
precessional effects. We are interested here in the rate of 
change of J, and therefore select the projection of J a in 
the direction of s a . This is 



J 



e abcQ a p£ Pb S ° 



(4.8) 



When we insert Eq. (3.29) in this we find once more 
that the Q? H C ab and £ ab terms make no contribution. We 
obtain 

j = -^eVM4s + Ma) ■ (4-9) 



We relate this to the rate of change of the body's rota- 
tional energy E rot = \lVt 2 H: where I is the body's mo- 
ment of inertia. With J — I fin we find that E mt = J. 



Taking into account Eq. (4.4 1, Eq. (4.9 1 becomes 

^A£*\£- al + A£- al ). (4.10) 



E IC 



n H j - 



Again we find that the effect vanishes in the absence of 
viscosity: It is the misalignment between the directions 
of the tidal forces and tidal bulge that allows the forces 
to exert a torque on the body. 

For an external object of mass M' on a circular orbit 
of radius b and angular velocity f2, the rate of change of 
the body's rotational energy is 



a 5 r [GM') 
18G W 



l\2 



n H (n-n H ). (4.ii) 



There is no change if the external object is corotating 
with the body. 



C. Heat dissipation 

The rate at which heat is dissipated within the fluid is 
given by (see, for example, Sec. 16 of Ref. (15"] ) 

Q=\pv J (d- a v h + d- b v & ) (d a v h + dV) d 3 x. (4.12) 
The velocity field in the norotating frame is obtained by 



transforming Eq. (3.201 using the rules spelled out in 
Sec. Ill E. We have 



t H e - b£ s x 

•2(4: + 



4GM 

■'be. 



'-)(£% + A£^ 



A£- h ,)x a x b x 5 



(4.13) 
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and this can be inserted within our expression for Q. 
A straightforward computation returns 

Q=^(£ al + A£ ah )(£- ah + A£- ah ). (4.14) 

Notice that this is positive-definite, as it should be. 

For an external object of mass M' on a circular orbit 
of radius b and angular velocity f2, the rate at which heat 
is dissipated within the fluid is 



and its averaged kinetic viscosity v. For an incompress- 
ible fluid with a uniform viscosity n 2 = &2 = § and 
t = ^-av / (GM) . [In the generalization I chose to pro- 
mote t to a dimensionful "Love quantity." An alternative 
would have been to define a new Love number v 2 and to 
generalize r as v 2 avj{GM). The former strategy is more 
useful for the purpose of establishing a correspondence 
with the black-hole results.] Similarly the generalization 
of Eq. (13.29} appears in Eq. O I. 



• _ qBt {GM'f 



D. Energy balance 



The expressions obtained for W and J — see Eqs. (4.3 I 
and (4.10), respectively follow directly from the ba- 

and the relation of 



( 4 - 15 ) sic definitions of Eqs. dOf, C 



Eq. (3.29 I between Q ab and £ a j,. The definitions continue 



to apply even after Eq. (3.29 I is generalized to Eq. ( 1.6 1, 
and we arrive at 



The work done by the tidal forces is used partially to 
increase the body's rotational energy, and is partially dis- 
sipated into heat by the fluid's viscosity. Energy balance 
requires that 



W 



(fc 2 r)a 5 
3G 



£ ab {£- al + A£- al ) 



and 



W = Q + fl H J, 



(4.16) 



n H j = 



(k 2 T)a l 
3G 



A£ al {£- a i 



A£- ai ) 



(4.17) 



(4.18) 



and we see from Eqs. (4.3 1, (4.101, and (4.141 that the These are Eqs. (1.8 1 and (1.9 1, respectively (in which 



equation is indeed satisfied. 



E. Generalization to arbitrary internal structure 
and correspondence with black-hole physics 

The results derived previously apply specifically to an 
incompressible fluid. It is easy, however, to generalize 
them so that they apply to a fluid body with an arbitrary 
internal structure (generated by an arbitrary equation of 
state) . The trick is to preserve the main structure of each 
equation, including the scaling with each dimensionful 
quantity, but to insert various dimensionless coefficients 
where appropriate. These coefficients are called "Love 
numbers," and their role is precisely to incorporate the 
details of the body's internal structure. This strategy is 
viable so long as one is interested only in the structure 
and scalings of the main equations, and one is satisfied 
not to know the exact value of each Love number. A 
detailed computation of the Love numbers for a selected 
equation of state would require the hard work of repeat- 
ing the analysis of fluid perturbations for this equation 
of state. 

The first result that can be generalized in this way is 
Eq. (3.21), which relates Qjk, the quadrupole moment 
of the body's mass distribution, to £jk, the quadrupole 
moment of the tidal field. The generalization appears in 
Eq. ( 1.5 ); it involves a rotational Love number n 2 , a tidal 
Love number fc 2 , and a viscous delay r that is known to 



the overbar is removed from the indices to simplify the 
notation, and G is set equal to unity). We can then 



invoke the energy-balance equation (4.161 to motivate 
the generalized form of Eq. (4.141, 



Q = {k2T)a5 (£-«» - 
* 3G V 



A£ ab )(£ ah + A£- a - b ) 



(4.19) 



Equations (4.171, (4.18), and (4.19) can be compared 



scale as av/(GM) with the body radius a, its mass M, 



with the corresponding black-hole results displayed in 
Sec. I. We associate W with M, J with J, and Q 
with (k/8it)A, where k is the black hole's surface grav- 
ity and A its surface area. Wit h thi s association the 
energy-balance equation of Eq. (4.161 becomes M — 
(k/8tt)A + Uff J, the first law of black-hole mechanics. 
The details of this comparison were already described in 
Sec. I. 
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